X-643-67-510 


M 


lid X- 35 9 74 


ITERATIVE INTERPOLATION AND 
NUMERICAL DIFFERENTIATION 
USING DIVIDED DIFFERENCES 


E. R. LANCASTER 


OCTOBER 1967 



NASA 




GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 

67-39691 


(ACCESSION NUMBER) 

z 

<THR^) 

(PAGES) . 

(NASA CR OR TMX OR AD NUMBER) 

_~7f 



ITERATIVE INTERPOLATION AND NUMERICAL 
DIFFERENTIATION USING DIVIDED DIFFERENCES 


E. R. Lancaster 


October 1967 

Goddard Space Flight Center 
Greenbelt , Maryland 




ITERATIVE INTERPOLATION AND NUMERICAL 


DIFFERENTIATION USING DIVIDED DIFFERENCES 
E. R. Lancaster 
ABSTRACT 

An iterative method, of interpolation based on divided differences 
is developed which is about three times as fast as the iterative methods 
of Aitken and Neville. The method is extended to numerical differentiation, 
with similar improvements over Hunter's extension of the Aitken and 
Neville methods. 
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ITERATIVE INTERPOLATION AND NUMERICAL 
DIFFERENTIATION USING DIVIDED DIFFERENCES 

Given values of a function f for Xq, Xj., x 2 divided dif- 

ferences are computed from 




d. = f( 
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x i } ; 
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We can 

arrange 

these values 

in a table 

as follows: 

*o 
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iH 
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x 2 

$-20 

&21 

d-22 



*3 

^30 

^31 

^■32 

^-33 


x 4 

^40 

d 41 

^42 

d-43 

d 44 


Letting d. = d , the polynomial which takes on the values 
J J J 

f(x ) at x for k = 0, 1, 2, . . . , i is given by 
k K 
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i 

V x) = I d j Pj <*), 

J=0 

P Q (x) a 1, 

Pj(x) = (x-Xo) (x-xi)...(x-x J ;L )i J- 1, 2, ..., i. 

This can be written 

y i (x) " y i-i (x ) + d i Pi (x )» i - 1, 2, 3, ... 

P i (x) a (x-x i _ 1 )p ii rx), 1-1, 2, 3, ... 
y 0 (x) a do, P Q (x) a 1 

This is a convenient iterative method for interpolating f at x 
It is approximately three times as fast as the methods of Aitken [l] 
and Neville [2] . 

We rewrite y ± (x) and p (x) in the forms 

i 

Pl (x) . I b 4J x J , 

J“0 
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y i (x) - l c u 

J=0 
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and tabulate the coefficients as follows: 
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b 2 2 
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C 31 

C 32 
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c 40 

C 41 

c 42 

C 43 
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i = 0, 1, 2, 


c . . 

1J 


c . _ . + d. b. . : 

i-i, j i ij" 


i = 1, 2, 3, 


b 44 


2, 3, . • • ) j — 1, 2, • . • , i-1 
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y j 0 j • • • y i “ 1 
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The^ith row of the c -table gives the coefficients of the poly- 


nomial 


2 with values f(x^) at x * x fc ; k = 0, 1, . .., i. 


j=0 

If the origin is placed at the point of interpolation, then c iQ is 
the interpolated value based on the points (x^ f(x^)); k = 0, 1, . .., i. 
m.' c^ m is an approximation of the mth derivative at x = 0, i.e., 
f (m) ( 0 ), based on the same points. If approximations of f^(o) are 
not desired for j > m, then it is not necessary to compute b. . and 

1 J 

c^ for j > m. 

This method of iterative numerical differentiation is considerably 
faster than that advocated by Hunter [3], which is similar to the 
methods of Aitken and Neville for iterative interpolation. 

The above algorithm can also be used to solve a set of Vandermonde 
equations or to invert a Vandermonde matrix. 

The iterated divided difference method of this paper was programmed 
by C. R. Herron and applied to the example given by Hunter [3]. The 
results were identical. 
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